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Abstract 



Cohesive zone models provide an illuminating and tractable way to include constitutive non- 
linearity into continuum models of defects. Powerful insights have been gained by studying 
both dislocations and cracks using such analyses. Recent work has shown that as a result of 
the locality assumption present in such cohesive zone models, significant errors can be made 
in the treatment of defect energies. This paper aims to construct a non-local version of the 
Peierls-Nabarro model in which the atomic level stresses induced at the slip plane depend in 
a non-local way on the slip degrees of freedom. The non-local interplanar kernel is computed 
directly from atomistics and is used to evaluate both the structure and energetics of planar 
dislocations. The non-local formulation does not significantly change the dislocation core 
structure from that obtained with the local model, but the new formulation leads to signif- 
icant improvements in the description of dislocation energetics for dislocations with planar 
cores. 

1 Introduction 

The appropriate method to use in modeling the behaviour of a material is often dictated 
by some intrinsic length scale in the problem. When considering atomic scale defects, the 
nanometer length scale is dominant, and fully atomistic models are often required. On 
the other hand, macroscopic models are concerned with bulk properties of specimens with 
dimensions on the micron scale or larger. In this regime, the simplifying assumptions of 
continuum mechanics are justified. Recently, there has been increased interest in modeling 
the so-called mesoscale regime, the range of length scales that are often too large for fully 
atomistic models but too small for discrete lattice effects to be ignored. In this regime, 
neither atomistic modeling nor continuum mechanics is satisfactory, and new models which 
incorporate features from both approaches seem to be necessary. 

One class of model which serves as a bridge between the microscopic and macroscopic 
approaches is that of cohesive zone models. Using a cohesive zone model allows bulk regions 
to be handled using conventional continuum mechanics, while atomistic effects are incorpo- 
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rated only at certain interfaces where it is deemed that they are important. Models of this 
type were first proposed by Peierls to describe dislocations, by Barenblatt to model 
fracture processes and by Dugdale [f§, and Bilby, Cottrell and Swinden || to estimate plastic 
zone sizes ahead of cracks. A review of a number of examples of this approach may be found 
in Miller and Phillips |J, while this paper focuses on the specific example of the Peierls 
dislocation. 

The Peierls dislocation model has received renewed interest of late, being used, for ex- 
ample, by Rice || in his description of the brittle versus ductile behaviour of crystals. The 
attractiveness of the Peierls framework is that it offers an analytically tractable (or at least 
numerically expedient) continuum model which incorporates nonlinear features resulting 
from the presence of the discrete lattice. The model assumes that the atomistic features of 
a dislocation are confined to a single atomic plane referred to as the slip plane. It is only at 
the slip plane that discontinuities in the displacement fields are assumed to take place, and a 
special constitutive law is used to account for them. Away from the slip plane, the remainder 
of the bulk crystal is treated as a linear elastic medium. In conventional treatments, the 
constitutive law at the slip plane is simplified by assuming that the energy and stress depend 
only locally on the slip distribution, despite the non-local nature of atomic interactions. In 
effect, while the Peierls model does incorporate a periodic length into the local relation be- 
tween stress and slip, it does not define a length scale below which the non-local effects due 
to slip gradients become important. 

In earlier work || it has been demonstrated that the failure to incorporate non-local 
effects is a serious concern, since full atomistic calculations show that the gradients in the 
slip distribution that are present in realistic dislocation cores are often too large for non-local 
effects to be ignored. By postulating a new form for the energy of a slip distribution, we 
formulate a model which includes a non-local term while at the same time reduces to the 
original Peierls model in the limit where slip gradients become small. 

Many examples exist where a certain model, derived under specific assumptions, is pushed 
to the extremes of its range of applicability and as a result must be corrected through 
inclusion of higher order terms. Anharmonic models of lattice vibrations |7[], ||, recent 
work on gradient models @,[10], and several papers by Eringen and co-workers [fLlJ|- ||13|| on 



non-local continuum theories are but a few examples of this strategy. Similar extensions 
are made here in the context of the Peierls framework. Due to the inclusion of a new term 
in the expression for the energy of a slip distribution, the model is now able to capture 
non-local effects due to the presence of gradients in the slip distribution. This improvement 
is demonstrated via comparisons to fully atomistic calculations of slip plane energies, and is 
then used to model dislocation core structures. 

In section 2 we begin with a brief description of the original Peierls framework and its 
breakdown in the limit of rapidly varying slip distributions. In section 3, we present the 
details of our non-local formulation and describe how atomistic calculations are used to 
build the necessary non-local constitutive model. We then use the non-local Peierls model 
to compute the energetics of crystalline slip, and show that the non-local model is in better 
agreement with purely atomistic results than was the classical Peierls framework. Finally, in 
section 4, the non-local model is used to obtain the core structure of a (100) [Oil] dislocation 
in fee Al. 



2 Breakdown of the Traditional Peierls Framework 

In this section, we present a brief explanation of the Peierls model, but focus mainly on 
describing the breakdown of the local cohesive zone assumptions. For a complete description 



of the Peierls model, see Hirth and Lothe 14 



The Peierls model assumes that a dislocation can be described as two elastic half-spaces 
joined at a common plane on which there is a discontinuous jump in the displacement fields. 
We adopt the convention that the slip plane is the x — y plane, with the dislocation line 
along the y-axis. The discontinuity in displacements due to the presence of the dislocation 
is referred to as the slip distribution, d(x,y) = u + (x,y) — u~(x, y), where ir^a;, y) are the 
displacement fields just above and below the slip plane. We confine our discussion to the 
simple case of plane strain in the x — z plane, and to the situation where only one component 
of 5 is non-zero. This allows us to write 5 as a scalar, 5(x). Dealing with the more general 
case is a straightforward extension of the results presented below, while the simple case being 
discussed here is better suited to demonstrate our arguments. It is important to bear in mind 
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that S(x) is not constant as in rigid slip, nor is it a simple step function as in the Volterra 
model of a dislocation. Rather, the slip distribution varies from zero at a point on the slip 
plane far from the dislocation core to a full Burgers vector once the core is traversed. The 
slip distribution is assumed to lead to atomic level forces due to the interaction between the 
slipped surfaces, thus providing the tractions, t(x), on the elastic regions. These tractions 
take a simple form that depends only on the local slip discontinuity 

t(x)=t[5(x)], (1) 

and can be determined on the basis of atomistic calculations. Early work within the Peierls 
framework assumed a simple periodic form for eqn. (fj), with the periodicity tied to the 
Burgers vector (see, for example, Foreman, Jaswon and Wood ||15|| ). More recently, highly 
accurate atomistic calculations have allowed for the direct calculation of t(8), improving 
the agreement between models built upon the Peierls framework and the results of direct 



atomistic simulation (see, for example, Rice, Beltz and Sun |L6] ). 

For the purposes of this paper, we find it more convenient to consider dislocation ener- 
getics rather than the resulting forces and tractions. According to the Peierls model, the 
energy of a dislocation is made up of two parts - the elastic energy contained in bulk regions 
and the misfit energy associated with the slip plane. The first component, the elastic energy, 
is fully defined once the elastic constitutive law for the bulk regions is specified, while the 
misfit energy is computed as an integral over the slip plane. While this is generally a surface 
integral, for the simplified geometry of a straight dislocation oriented along the y-axis, the 
integral reduces to 



E L = / $[6(x)]dx, (2) 



-oo 



where is understood to be the local misfit energy per unit length along the dislocation 
line. This convention will be maintained throughout the paper. In eqn. (|2[), $(£) is given by 

s 

$(<5) = J T (5')d5'. (3) 
o 

$(<5) is referred to as the interplanar slip potential, and can be thought of as the energy 
cost associated with slipping one block of atoms over another by an amount 5. As with the 



tractions t(8), the interplanar slip potential can also be obtained through simple atomistic 
calculations (see, for example, Rice et al. [16| or Kaxiras and Duesbery jT^I ). 



The energetic description provided above introduces an important assumption that has 
been central to the Peierls framework and results in a formulation that is strictly local. 
Specifically, this assumption arises from the fact that the interplanar potential is computed 
on the basis of a purely uniform slip distribution, even though the actual slip distributions 
of interest are non-uniform. The energy of the non-uniform slip distribution is found using 
eqn. (0), which effectively divides the slip distribution into infinitesimal slip steps, samples 
the interplanar potential for each of these steps, and sums the results. The tacit assumption 
of this approach is that despite the non-uniform nature of the slip distribution, the local 
environment at each point can be considered to be approximately uniform, allowing for the 
local slip energy to be determined from the interplanar slip potential. It is expected that 
as long as the gradients in the slip distribution are small, this approximation will be valid, 
but as the gradients become more severe, the approximation will breakdown. One result 
of interest is to quantify what the maximum acceptable gradients are, and to determine 
whether or not the slip gradients occurring in real dislocation core structures exceed these 
values. 

These questions about the validity of the assumptions in eqn. (|2]) were addressed by Miller 
and Phillips 0. The breakdown of the local assumptions was quantified by computing 
the energy of a number of idealized slip distributions in two ways, first using atomistics 
(considered "exact" in this context) and then using the approximation embodied in eqn. @. 
The parameters in the slip distribution were then selectively varied, allowing for control of 
the severity of the slip gradients and a direct comparison of the two methods of obtaining 
the slip plane energy. These calculations demonstrated a clear breakdown in the locality 
assumption. Further, it was found that the gradients associated with the slip distributions 
for simulated core structures are of the same order of magnitude as those for which eqn. (Q) 
failed. 

In fig. (|^), we reproduce the results of Miller and Phillips for the {001}[110] slip system 
in fee Al and include new results for the {111} [110] slip system. The figure shows a plot 
of misfit energy per unit area of the slip plane for various slip distributions. For periodic 
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distributions, this is computed by dividing the misfit energy per period by the periodic 
length. For slip distributions which are not periodic, an effective slip area was used, which 
is defined as the area over which the slip is greater than 1% of its maximum value. The 
following slip distributions were used, 

5(x) = Asin^-, (4) 
w 

6(x) = AeM-\0% ^ 
For these slip distributions, w can be varied in order to control the slip gradients. Small w 
corresponds to large gradients, and hence the regime in which we expect the local approx- 
imation to fail. To build these slip distributions for an atomistic calculation, one divides a 
crystal in two, and imposes a different deformation field in the upper and lower halves. Ju- 
dicious choice of these deformations leads to the appropriate discontinuity at their common 
surface, and allows for the computation of the misfit energy due to that discontinuity. This 
is the exact atomistic energy which is compared to the results of eqn. (0) in fig. ([[]). Further 
detail of the atomistic calculation of misfit energy is given in section 3. 

Figure (0) includes two different types of slip deformation. In the slip distributions asso- 
ciated with straight dislocations, two types of gradient effects are possible. First, both the 
direction in which the slip distribution is changing and the direction of the slip discontinuity 
itself can be the same. This is characteristic of the slip distributions for pure edge dislo- 
cations, and therefore we refer to any slip distribution for which the 5 vector and the V5 
vector are parallel as an "edge" -type distribution. The second type of gradient effect occurs 
when the direction of the slip discontinuity is perpendicular to the direction along which it 
is changing. This is characteristic of a pure screw dislocation, and therefore we refer to such 
an instance as a "screw" -type slip distribution. It is possible to mix these two effects, but 
for our present purposes we consider only the pure edge and pure screw cases. These two 
types of slip distributions are analogous to longitudinal and transverse phonons. 

The results in fig. (0) demonstrate the breakdown of the Peierls assumption for the slip 
distributions of eqns. (D and (D for both the {001}[110] edge distribution and the {111}[110] 
screw distribution. The plots give the energy of the slip distributions as a function of the 
parameter w. In all of the plots, the prediction of the local slip approximation for the 
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misfit energy is shown as the constant dashed line near the top of each graph, a striking 
manifestation of the lack of a characteristic length scale in the local cohesive zone approach. 
For periodic slip distributions, it is easy to see why the Peierls model approximation is 
independent of the parameter w. The plots of fig. ([[]) are of the quantity 



££ = 1 

W W 



J $[v4sin(27r x/w))dx, (6) 



where E\ is the energy of a single periodic length of the slip plane. The fact that this 
integral correctly represents the energy of a single periodic length is discussed in section 3. 
By making the change of variables y = x/w, one can see that the expression is independent 
of w. The data points in fig. (|1|) come from atomistic calculations of the same misfit energies. 
Note the divergence from the local estimate at small values of w. 

The first mission of the non-local model will be to improve the agreement between the 
atomistic and continuum predictions for the slip energy shown in fig. ([I]). The solid line in 
this figure shows the improvements made by using the non-local model, and will be discussed 
in more detail at the end of section 3. 



3 The Non- Local Formulation 

The notion that any field variable (for example stress, strain or temperature) is known 
pointwise, and depends only on other state variables at that point is a natural conclusion of 
observing the macroscopic behaviour of materials. The result is the assumption of locality - 
one of the basic tenets of the Peierls model, most other cohesive zone models, and classical 
continuum mechanics in general. However, on the atomic scale the state of each atom is 
influenced by a finite cluster of its neighbours, not only by the state at an individual atomic 
site. The non-local formulation of the Peierls model outlined in this section is proposed as a 
simple way to include the non-local nature of atomic interactions into cohesive zone models 
such as the Peierls framework. 

The modification made to the Peierls model in going from the local to non-local formu- 
lation involves the inclusion of a non-local contribution to the misfit energy. The addition 
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Figure 1: Energy of idealized slip distributions as a function of the parameter w, which 
determines the gradients in the slip distribution. These results demonstrate the failure of 
the local slip approximation for small w. Exact atomistic energy (filled circles); local model 
(dashed line); non-local model (solid line). 
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of this term has far-reaching implications, though it is the only conceptual change that is 
introduced. 

Consider again the misfit energy of a dislocation as given in the local framework of 
eqn. (0). Non-locality is incorporated by the addition of a term which should vanish for 
slowly varying slip distributions while at the same time capturing gradient effects when they 
are present. For a straight dislocation, we postulate the non-local misfit energy to be given 
by 

oo oo oo 

E MF = J $[5(x)}dx + J J K(x - x')5(x)5(x')dx dx'. (7) 

— oo — oo — oo 

The additional term contains the non-local interplanar kernel K(x — x'), which weights the 
non-local contributions to the total energy. On physical grounds, we assume that K(x—x f ) = 
K(x' — x), or that the influence that one point has on another depends only on the distance 
between these points. The new model requires only that we determine a suitable form for the 
non-local kernel K(x — x'). We propose to adopt a similar strategy to that used in obtaining 
$(5), whereby we extract a numerical reckoning of K(x — x') from atomistic calculations and 
then fit these results to simple analytic forms. In the past, non-local expressions of this type 
have been advanced without the benefit of atomistic calculations to determine the influence 
function It seems possible that the methods presented here can be adapted to those 



cases as well. 

Given the non-local kernel, the traction at the slip plane can be written as the first 
variation of the misfit energy functional with respect to the slip distribution, yielding 

dE MF d<f>[8(x)' 



05 35 



2 / K(x - x')5(x') dx'. 



-oo 



Thus, for the non-local model, the traction r at any point on the slip plane depends on the 
entire slip distribution, whereas in the local formulation the stress is determined pointwise. 
It is this additional feature of the new model which allows it to capture the gradient effects 
discussed earlier. 

To determine the non-local kernel, we solve eqn. (0) for K(x — x') by imposing a number 
of special slip distributions, and use the non-local term to exactly fit the energies of these 
slip distributions to the atomistic result. This procedure amounts to the assumption that 
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K(x—x r ) is independent of the form of the slip distribution. While this may not be rigorously 
true, it is postulated that the non-local kernel obtained in this way will be reasonably effective 
for general slip distributions. Based on this assumption, it is convenient to choose slip 
distributions which allow for a determination of the Fourier components of K(x — x'). For 
the purposes of our calculation we place two demands on the slip distributions. First, that 
it be periodic with period L, and second, that the maximum amplitude of the distribution 
is small with respect to the lattice parameter of our material. The first requirement is one 
of computational convenience, allowing for the use of periodic boundary conditions when 
determining the exact atomistic energy of a slip distribution, while the second requirement 
is made in order to allow us to rewrite the interplanar potential, $, using its quadratic 
approximation. This step proves essential in determining Kix — x') in much the same way 
that one uses the quadratic approximation to match the shear modulus in the case of the 
Frenkel sinusoid model. It simplifies the local term in the energy for the purposes of finding 
the non-local kernel, but once the kernel is computed the exact form of the local term will 
be reinstated. This ensures that the non-local model will still be applicable to large slip 
deformations. 

The energy expression of eqn. (|?p yields the energy of the slip distribution over the entire 
slip plane. On the other hand, atomistic models naturally provide us with a way to compute 
the energy of only a finite section of the slip plane. If periodic boundary conditions are 
used in the atomistic model, as will be used here, then the energy obtained from atomistics 
is the energy of a single periodic length of the slip distribution. In order to make valid 
comparisons between the exact atomistic misfit energy and the energy obtained from the 
non-local model, we must find the non-local expression for the energy of a single periodic 
length of slip. Equation (!?]) can be re-written in the form 



Emf(5)= £ 



(n+l)£+£ (n+l)L+$ oo 

J c5(x) 2 dx + J J K{x — x')5(x)5(x')dx' dx 



(9) 



where the constant c arises from the treatment of the interplanar potential via its quadratic 
approximation and £ is some origin where < £ < L. It is easy to show that for periodic 
slip distributions, the integrals inside this sum are independent of n and £, and hence each 
periodic length in the variable x contributes the same amount to the total misfit energy. We 
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can then write the misfit energy of a single period of the slip plane, E® IF (5), as 



L oo 



E° MF {8) = Jc5(x) 2 dx + J J K(x - x')S(x)5(x')dx' dx, (10) 



-oo 

,-0 



and note that for a given slip distribution, E MF {5) is a quantity which can be computed 
directly from atomistics. For this purpose, it is convenient to take S(x) to be 

5(x) = S q (x) = Asmqx, (11) 

where q = 27r/L and A is much smaller than the lattice constant of the crystal under 
consideration. We insert this form of S(x) into the expression for E° MF . Making a change of 
variables z = x — x' and recalling that K(z) = K(—z), the energy expression becomes 

cA 2 L „ r L 

dz dx, (12) 



cA Ij f n r°° 

E° MF (6 q ) = — + A 2 singx / K(z 

l JO J-oo 



2i 

where we have made use of the exponential form of the sine function. We recall the definition 
of the Fourier transform, f(k), of a function f(x) as 

oo 

/(A) = / f{x)e- ik *dx, (13) 

— oo 

and also the inverse Fourier transform 

oo 

/(*) =7T I fWe^dk, (14) 

— oo 

and use eqn. flT3"|) to write the non-local energy expression in terms of the Fourier transform 

of the non-local kernel, K(q). Noting again that K(z) = K(—z), we see that K(q) = K(—q), 

and the expression for the total energy becomes 

rA 2 T r L 
E° MF (S q ) = + A 2 K(q) / sin 2 qxdx. (15) 

2 Jo 

Evaluating the integral leads to an explicit expression for the non-local interplanar kernel in 
Fourier space; 

K{q) = %P - c. (16) 

At this point, the strategy is to obtain the energy dependence on the Fourier variable q 
numerically through an atomistic model. 
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3.1 Atomistic Determination of the Non-Local Interplanar Kernel 

The procedure for determining the energy of a given slip distribution was outlined in section 
2 and is described in detail by Miller and Phillips ||. Here, we use the same computational 
approach which is described briefly below. 

We desire the energy of the slip distribution given by eqn. ([11]) as a function of the pa- 
rameter q. To make the computation, the crystal is divided into an upper and a lower half. 
By imposing the appropriate displacement fields on the two half spaces, we can create the 
desired slip distribution at their common interface. For example, consider the slip distribu- 
tion of eqn. (|Tl|). For the case of edge-type slip, the slip distribution can be obtained by 
imposing the following displacement fields on the crystal: 

u + = [(A/2) sin gx, 0,0], (17) 

u = [(-A/2) sings, 0,0], (18) 

where u + represents the displacement field in the upper half of the crystal and u~ is that 
in the lower half. This combination of displacement fields leads to a slip distribution of the 
form 

5(x) = — u~ = Asmqx (19) 

at the plane where the upper and lower half spaces meet. 

Using any convenient atomistic model (we have used the embedded atom method (EAM), 
see |L8| for example.) it is then possible to compute the energy of this crystal. The result of 
such an atomistic calculation is an energy consisting of two parts, that due to the interface 
and an additional elastic strain energy due to the deformation in the bulk regions of the 
crystal. The elastic energy can be found directly by computing the energy of each half space 
separately, using a periodic computational cell containing no slip discontinuities. This energy 
can then be subtracted from the total energy of the configuration which includes the slip 
jump, leaving the misfit energy, E^ F . 

By repeating this procedure for a sequence of values of q, we obtain a discrete represen- 
tation of the non-local kernel K(q) in Fourier space. Results of such calculations using the 
EAM potentials for Al of Ercolessi and Adams [^] are presented in fig. (^). The various slip 
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systems are either edge, screw or "mixed", where edge and screw are as defined in section 
2. The "mixed" system is in this case the direction associated with the 30 degree Shockley 
partial. All of the curves are given in non-dimensional form using the constants given in table 
([[]), where 7 MS is the unstable stacking fault energy for a given slip system. The values in 
the table are obtained by rigidly sliding two blocks of atoms with respect to one another and 
allowing relaxations in only the out-of-plane direction. The remaining constants in the table 
are b, the Burgers vector, and d, which is defined below. The representative examples in 
fig. (0) demonstrate a number of the characteristics of the non-local interplanar kernel when 
constructed in Fourier space. All such functions are even and periodic, and therefore only 
the first half period of each is shown. The periodic length of K depends on the slip system 
being considered, and is given by 2n /d where d is the distance between planes with normal 
parallel to the direction of the slip gradient vector. Note that these planes are a set of y-z 
planes, which are perpendicular to the x-y slip plane. Another characteristic feature of the K 
function is that it is exactly zero at q = 0, which implies that in the limit of very slowly vary- 
ing slip distributions, the non-local correction vanishes and we recover the classical Peierls 
framework. It is also worth noting that while most of the K functions are negative, we see 
that for the case of screw-type slip distributions on the {001} [010] slip system K is entirely 
positive. This suggests that there is no set rule about whether the non-local contribution 
to the misfit energy is positive or negative - it will depend on both the slip system being 
considered and on the Fourier components of the slip distribution itself. Comparing the var- 
ious curves in this figure can help us to understand the relative importance of the non-local 
effect for various slip systems. Note that the strongest effect is associated with the (100) [010] 
edge system. It is interesting that this same slip system, but in the screw orientation, shows 
the weakest non-local correction, and that the sign of this correction for the various Fourier 
components of the slip distribution is reversed. This is because the screw-type sinusoidal 
slip distribution represents a more severe misfit configuration than the edge-type for this 
slip system. Meanwhile, the local model predicts the same energy for either edge or screw 
slip due to the symmetries of the {100} planes. For this slip system, the local prediction for 
misfit energy somewhat overestimates the edge energy, while somewhat underestimating its 
screw counterpart. The result is this seemingly anomalous behaviour of the screw-type slip. 
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Figure 2: Examples of the Fourier space representation of the non-local interplanar kernel, 
K, for various slip systems in Al, as obtained using the EAM potentials of Ercolessi and 
Adams [19]. See table (1) for normalization constants. 



Slip System 


b(k) 


d/b 


lus (eV/A 2 ) 


(100) [010] screw 


4.032 


1/2 


0.04065 


(111) [110] screw 


2.851 




0.02705 


(111) [211] mixed 


1.646 




0.008032 


(100) [011] edge 


2.851 


1/2 


0.02521 


(100) [010] edge 


4.032 


1/2 


0.04065 



Table 1: Burgers vector b, interplanar spacing d and unstable stacking fault energy 7 US for 
the slip systems considered in fig. (2) as obtained from the EAM Al potentials of Ercolessi 
and Adams [19]. 
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3.2 Approximate Analytic Form in Real Space 

Numerical inversion of the Fourier representation of K deduced above can be problematic 
due to the fact that K is known only for a discrete set of points in Fourier space. Therefore, 
we propose to fit K in Fourier space with a cosine series, which can be easily transformed 
into real space in the form of a sum of Dirac delta functions. Because K(x — x') exists 
as part of an integrand over the entire slip plane, these Dirac deltas will have the effect of 
reducing the dimension of the integral in eqn. (|7|). We have found that excellent fits of the 
numerical data can be obtained using only the first few terms of such a cosine series, and 
one can therefore write 

a N 

K{q) = + cos ndg, (20) 

1 n=l 

where a n are the fitting parameters and d is the distance between planes perpendicular to 
the slip plane, as described previously. The number of fitting parameters N will depend on 
the range of the atomistic potentials used in obtaining K, but we have found that for the 
potentials and slip systems considered here, the value of N typically does not need to exceed 
five. This form for the non-local kernel allows for easy Fourier inversion, yielding 

a N a 

K(x — x') = -^-Sd(x — x') + ^2 "^"(^-d( x — x' + nd) + $d( x — x' — nd)\ (21) 
2 n =i ^ 

where 3d is the Dirac delta function. This expression for the non-local kernel can be used 
in the original definition for the misfit energy, as well as in the expression for the tractions 
on the slip plane. The results are 

E MF = J •[*(,)]*, - £ / aM*Y - + (22 ) 

-oo n=1 -oo 

T (x) = [ "l n - V a n {25{x) - 6(x + nd) - 6(x - nd)), (23) 

96 n=l 

where we have eliminated from the expressions by using the fact that K{0) = 0. The 
S appearing in this expression is the slip distribution, and not the Dirac delta function Sd- 
The computed values of a n are plotted in fig. (|]) for the five slip systems considered here. 
In this figure, increasing n corresponds to sampling the slip distribution farther and farther 
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Figure 3: Non-local influence coefficients, a n , for the slip systems considered in fig. (2). 



from the point at which the energy and stress is being computed. It is clear that the non- 
local effects decays quite rapidly, and in all cases a\ (corresponding to the near neighbour 
non-local influence) is by far the most important contribution. Non-local influence is felt as 
far as five neighbours away in some systems, but beyond 05, the coefficients are negligible. 

Note that the expressions for the energy and stress are consistent with our intuition 
as to the nature of the non-local terms. First, in the limit when 5(x) is slowly varying, 
5(x ± nd) — > 5(x) and the above expressions reduce to the conventional local formulation. 
Second, the non-local effects for a given point arise as a result of sampling a discrete set 
of points along the rest of the slip distribution. These points are spaced by a distance d, 
which coincides with the spacing of atomic planes in the direction of the slip gradient vector. 
The discreteness of this sampling, together with the physical significance of the spacing 
of the sampling points results in a sensible atomistic extension of the original constitutive 
assumptions of the local model. 

Finally, it is interesting to note that the terms of the sum in eqn. can be viewed 
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as a linear combination of approximations to the derivatives of the slip distribution. The 
symmetry of our non-local kernel in Fourier space means that only even order derivatives 
contribute to this correction, and therefore it is possible to re-write eqn. (|23| ) in the form 

where the coefficients c n are the weights associated with the nth derivative, and 5^ is the 
central difference approximation to the nth derivative of the slip distribution. The c n are 
simply linear combinations of the a n already introduced. For example, if we take N = 2 in 
eqns. ( |2"B| ) and ([23]), and make use of the following central difference approximations for the 
second and fourth derivative of 8(x), 

d 2 5 8(x + d)- 25(x) + 8(x - d) 



dx 2 d 2 

d*5 _ 5(x + 2d) - A5(x + d) + 65(x) - A5(x - d) + 5(x - 2d) 
dx 1 ~ d 1 ' 

then we can solve for the coefficients c 2 and C4. In this case, these are found to be 



(25) 
(26) 



c 2 = -d 2 (a l + a 2 ), (27) 

c 4 = -d A a 2 . (28) 

This exercise highlights the parallels between our non-local approach and recent gradient 
correction models (e.g. [[|, |]10[), although in this paper we will not explore the gradient 
correction form of these equations further. 



3.3 Energy of Crystalline Slip Within the Non-Local Formulation. 

The simplest test of the non-local model is to re-examine the energetics of idealized slip 
distributions as originally presented in fig. (|J). Recall that in this figure the predictions of 
the local model are given by the constant dashed line near the top of each graph, whereas the 
data points are exact calculations of the misfit energy as determined using atomistics. The 
solid curves are the results of using the non-local model, which are in significantly better 
agreement with the exact results than are the local results. In each case, the energetics 



18 



using the non-local model exhibit the correct trend as a function of slip gradients, unlike the 
original model which does not capture gradient effects. The non-local model is not able to 
capture the anomalous upturn in the atomistic results for the {111} [110] screw orientation, 
a pathology of this particular slip deformation which results from putting atoms into highly 
unfavourable proximity across the slip plane. Nonetheless from the standpoint of purely 
energetic considerations, the non-local model shows a marked improvement over the local 
model, without any great cost in model complexity. 

Emboldened by our observations that the non-local correction improves the energetic 
description of crystalline slip, we now proceed to a result of greater interest, namely, the 
effect of the non-local terms in the context of realistic dislocation core structures. 

4 The Non-local Model of Realistic Dislocation Cores 

As an example, we apply the non-local model to the determination of the core structure for 
a straight dislocation with a planar core. For this purpose we consider the Lomer dislocation 
in fee Al, with a < 110 > line direction and Burgers vector | < 110 >. 

We obtain the Lomer core structure using two schemes. The first makes use of a simplified 
form for the interplanar potential $(5) that allows for a closed form analytic solution for 
the slip distribution in Fourier space. In the second determination of the Lomer core, the 
full-blown atomistic result for $(5) is used and the core structure is computed numerically. 
Each of these results is then compared to the slip distribution taken directly from the atomic 
positions resulting from full relaxation of the atomistic degrees of freedom. The first scheme, 
which admits of an analytic solution, is of interest as a method for testing the numerical 
procedures used in the second scheme. At the same time, it demonstrates that considerable 
analytical progress can be made with the non-local formulation. The second, fully numerical 
solution demonstrates that stable core structures are readily obtainable within the non-local 
formulation. We will see that while the effect of the non-local terms on the slip distribution 
and core structure are subtle to the eye, they represent a significant improvement when 
quantified in terms of the predicted misfit energy of the dislocation core. 
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4.1 Eigenstrains solution for 5 within the quadratic well approx- 
imation. 

In this section, we describe the procedure used in determining an analytic result for the 
core structure of the Lomer dislocation. The important approximation in this procedure in 
relation to the results of the subsequent section is the simplified form used for the interplanar 
potential $(5). 

Within the Peierls framework, the determination of 5(x) is equivalent to finding the core 
structure of the dislocation. Once 5(x) is found, the elastic displacement fields away from 
the slip plane can be computed from an integral of the Volterra kernel over the entire slip 
plane. As well, the misfit and elastic energy are then fully specified by the slip distribution. 
Therefore, we seek the slip distribution S(x) which minimizes the total energy functional. 
The energy associated with the dislocation, which is a functional of the slip distribution S(x), 
can be expressed as the sum of three parts 

Etot = El + Enl + Eb- (29) 

El is the misfit energy as determined within the traditional local model and given in eqn. @, 
Enl is the non-local correction term added to the misfit energy in eqn. (]?]), and Eb is the 
elastic energy of the bulk region. The elastic term is obtained by superimposing the elastic 
interaction energy for a distribution of infinitesimal dislocations with Burgers vector density 
—dd/dx, 

1 ,co ,co R d5(x)d5(x')j , 

EB= 2 oo oo El0g( b^ ) dx dx' dXdX - (30) 
Zj J —oo J — oo X — X LLX tlx 

In this equation, R is a measure of the size of the bulk region. Making use of integration by 
parts, it is possible to eliminate the constant R from the analysis. B is defined as 

B = 2C lj s i s j (31) 

where is the ith component of the slip direction for the dislocation of interest and Cy is 
the prelogarithmic energy tensor, discussed in detail, for example, in Bacon, Barnett and 



Scattergood |2(J. In the case of an isotropic solid B reduces to 



Bis ° = 2k(1-uY (32) 
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Figure 4: (a) Quadratic well approximation to the local interplanar potential, (b) Eigenslip, 
5 E (x) for the dipole configuration. 
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For the purposes of obtaining an analytic solution, we make the approximation that 
the local interplanar potential $(<5) takes the form of a periodic array of quadratic wells as 
shown in figure (f|a). We then employ the method of eigenstrains as described by Mura pi] 
to minimize the energy functional with respect to the slip distribution for a dislocation dipole 
with spacing r. The dipole configuration is used so that the slip distribution goes to zero at 
x = ±00. By making r large, the dipole becomes two isolated dislocations of opposite sign. 
This approach allows us to rewrite the local contribution to the energy as 

/oo 
c(5(x) - 8 E (x)) 2 dx. (33) 
-00 

In this equation the constant c is determined by fitting the quadratic wells to linear elasticity 
at small values of 8 — 8 E . This leads to c = /i/2a where a is the spacing between slip planes 
and \i is the relevant shear modulus for the slip system of interest. 8 E in this expression is 
the "eigenslip" which ensures that the dipole slip configuration is enforced during the energy 
minimization. 5 E (x) is shown in figure (|]b), where b is the Burgers vector and r is the dipole 
spacing. We will need the Fourier transform of this eigenslip, which is found to be 

r/2 9b kr 

b e - ikx dx = —sm—. (34) 

-r/2 k 2 V ^ 



Making use of eqn. (14]), we can replace 8, 8 and K with their Fourier transforms, allowing 



for a Fourier space representation of the total energy. Integrating by parts on the bulk energy 
terms and assuming that the order of integrations can be switched allows us to simplify the 
energy expression to 



tot 



$(*)*H0(f + ^ + f 1*1) - -Kk)S E (-k) + ^8 E (k)8 E (-k) 

Z7T 2tt 4 71 2tt 



dk. (35) 



Note that at this point, our solution is predicated upon a knowledge of the Fourier space 
features of the interplanar kernel. Taking variations with respect to 8 and setting 8E tot = 
we find 

S {k) = 2 ±^1 (36 ) 

where we have made use of the fact that K{k) is an even function and used eqn. ([53]) to 
replace the Fourier transform of the eigenslip. This expression gives us the Fourier transform 
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of the slip distribution for the dislocation dipole. Although it cannot be transformed into 
real space explicitly, the transformation integral is straightforward to perform numerically 
once the form of K is known. Thus we numerically compute the slip distribution from the 
expression 

s 2b r°° sin^fcosfcx „ , , 

6(x) = — / 1— dk. 37 

t Jo k[1 + mi + EM] 

This result is plotted in fig. ([5]) for the case of K = (which corresponds to the local Peierls 
model solution) and with K from fig. @ for the (100)[011]edge (Lomer) slip system. Note 
that only one member of the dipole pair is shown because of the symmetry of the distribution. 
We see in fig. (|5|a) that the effect of the non-local correction is to introduce small oscillations 
in the slip distribution and increase its slope slightly at the dislocation core. At first, the 
oscillations may appear to be unphysical, but recall that this continuous curve is really a 
representation of the slip for a discrete set of lattice sites. Therefore, only the values at the 
lattice sites are germane for the atomic positions implied by the solutions. It is interesting 
to note that this analytic model provides a reasonable approximation to the dislocation core, 
although the differences between the local and non-local results are hard to quantify. We 
should not expect the continuum model to accurately capture subtle details of the atomic 
core, but it is encouraging that the model does lead to a stable core configuration that is 
similar to the exact result. 

4.2 Numerical solution for the slip distribution. 

It is now of interest to see what effect the non-local correction has on the dislocation core 
structure when the correct atomistically obtained <&(5) is used. In this case, it is necessary 
to resort to a full numerical solution of the governing equation. 

The numerical approach for finding S (x) follows closely the approach of Beltz and Freund 



22fl , to which we refer the interested reader for more details. We begin with the well known 



Peierls-Nabarro equation |L4 



r[S(x)] = - , r d5{X ' ),d ? dx>, (31 



2tt(1 

into which we can substitute eqn. ( p3|) for r[5(x)]. By an appropriate change of variables, 
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Figure 5: Comparison of core structures for the Lomer dislocation in Al as obtained from the 
continuum model with exact atomistic results, (a) Structure obtained using the quadratic 
well model of the interplanar potential and (b) with the full numerical solution to the Peierls- 
Nabarro equation. 
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the domain of integration can be collapsed onto the finite domain (-1,1). This domain can be 
discretized and the integration carried out by making use of the Gauss-Chebyshev integration 
techniques described in Erdogan and Gupta f23fl . This approach reduces the problem to a 
set of nonlinear algebraic equations which can be solved iteratively via the Newton-Raphson 
method, leading to a discrete representation of the slip distribution 5(x). The results of such 
a calculation for the Lomer core structure are presented in fig. (]5|b). Again, we note that the 
resulting core structure is not changed dramatically relative to the local model as a result of 
adding the non-local correction. Neither the local nor the non-local model can fully capture 
the details of the core, and it is difficult, on the basis of visual inspection, to quantify the 
virtues of including the non-local correction when considering only the slip distributions. 

The fact that the structural differences between the local and non-local models are small 
does not imply that the non-local effect is unimportant. Instead, the important quantity to 
compare between the two models is their predictions for the total energy of the dislocation. 
The way the elastic energy is computed for the two models is identical, and the highly similar 
predicted cores mean that the value of the elastic energy for the two models will be about 
the same. On the other hand, the misfit energy is computed differently, and even identical 
core structures will lead to different misfit energies by virtue of the non-local energy term. 

Using the method described in Miller and Phillips f5[ to isolate the misfit energy from the 
strain energy in the bulk regions, we can compute the exact atomistic misfit energy of the 
lomer dislocation to be 0.1362eVyA 2 . This is accomplished by representing the atoms in the 
upper and lower bulk regions by nodes in a finite element mesh, and computing the strains 
(and consequent strain energy) in this mesh via the Cauchy-Born rule [23]. On the other 
hand, we can compute this energy using the local and non-local models from eqns. and 
(f2~2|) respectively. The local model predicts an energy of 0.1979eV/A 2 , about 45% greater 
than the exact energy. The non-local result is 0.1471el / /A 2 , only 8% greater than the 
atomistic result. This result demonstrates that although the core structures in the local and 
non-local formulations are for practical purposes indistinguishable, the non-local treatment 
of the core energies is significantly more accurate. We would expect similar adjustments to 
the cohesive zone model estimate of the Peierls stress. 
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5 Conclusions 



Motivated by previous work which demonstrated a failure of the local Peierls framework to 
accurately describe the energy of interplanar slip for the types of slip distributions found in 
real dislocation cores, we proposed a non-local formulation of this framework. We then pro- 
ceeded to outline a set of simple atomistic calculations whereby one can obtain the non-local 
interplanar kernel required in the formulation of the model. To demonstrate that the model 
improved estimates of the slip energy, we compared these results to purely atomistic calcu- 
lations and showed that the non-local model improves the agreement between the cohesive 
zone model energies and explicit atomistic energies. 

Given that the non-local model leads to better energetic descriptions of slip distributions, 
we proceeded to demonstrate the model by computing the structure of a straight dislocation 
with a planar core, namely, the Lomer dislocation in Al. It was found that the differences 
between the local and non-local results, and between either of them and the atomistic core, 
were subtle and difficult to quantify when considering only the spatial structure of the 
dislocation. However, the important quantitative measure of the performance of the models 
is the energy, which the local model predicted to be 45% larger than the exact atomistic 
result. The non-local model significantly improved the energetic description of the core, 
overestimating the exact result by only 8%. Nevertheless, the negligible changes in the non- 
local description of the core structure serve as a reminder that the cohesive zone approach 
appears to lack the flexibility to really serve as a generic basis for mixed atomistic and 
continuum studies of dislocations. Even in its non-local form, this framework restricts the 
slip to a particular slip plane, thus forbidding the emergence of complex cores such as those 
found in bcc metals. 

Future work in the context of cohesive zone models could include an attempt to integrate 
non-local effects into other models of the mechanical behaviour of materials. One example 
of interest is the model of Rice || and the numerous related works that describe dislocation 
emission from crack tips. Preliminary investigations in this area suggest that the concepts 
outlined here can be used for such problems . The basic equation to be solved in the case 
of a dislocation near a crack tip is a modified version of the Peierls-Nabarro equation with 
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additional terms due to the presence of the crack. Any realistic solution of this equation 
requires a numerical procedure, and the simple form of the non-local corrections in eqns. ( P2|) 



and (p3|) mean that the additional computational cost in the non-local formulation is small. 
The most serious obstacle in this case seems to be the proper treatment of the boundary 
conditions in the non-local setting, and most importantly the stress free boundaries at the 
crack faces. However, the fact that non-locality plays an important role at the atomic 
level means that a correct non-local treatment of the atomistically sharp crack may be an 
important contribution to our understanding of crack tip phenomena. 

Another possible direction for this work is the integration of the non-local formulation 



into the recently proposed "semi-discrete Peierls framework" of Bulatov [26]. Results using 
this version of the Peierls model show significant improvement to the conventional Peierls 
framework, but some error remains. It is possible that including non-local effects into this 
model may further improve its agreement with atomistics while retaining its tractability. 

A third future direction for this research should be an effort to compute other non-local 
kernels analogous to K(x — x'). For example, it may be possible to directly compute from 
atomistics the non-local elastic moduli introduced by Eringen and co-workers |H|— |[L3|| for 
bulk crystals. This would eliminate guesswork about their appropriate form and base them 



solidly on their atomistic underpinnings. Other recent work (P7||,[28|) substantiates our 
belief in the critical role played by constitutive non-locality in the description of atomic 
scale defects and calls for continued efforts to put such models on a clear analytic footing. 
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